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Abstract 

We study the classical and quantum evolution of a universe in which the matter source is a 
massive Dirac spinor field and the universe is described by a Bianchi type I metric. We focus 
attention on the those classical solutions that admit a degenerate metric in which the scale fac- 
tors have smooth behavior in transition from a Euclidean to a Lorentzian domain and show that 
this transition happens when the cosmological constant. A, is negative. The resulting quantum 
cosmology and the corresponding Wheeler-DeWitt equation are also studied and closed form ex- 
pressions for the wave functions of the universe is presented. We have shown that there is a close 
relationship between the quantum states and signature changing classical solutions, suggesting a 
mechanism for creation of a Lorentzian universe from a Euclidean region by a continuous change 
of signature. The quantum solutions also represent a quantization rule for the mass of the spinor 
field. 
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1 Introduction 



A question of great importance in cosmology is that of the initial conditions from which the universe 
has evolved. As is well known, standard cosmological models based on classical general relativity have 
no convincingly precise answer to this question. This can be traced to the fact that these models 
suffer from the presence of an initial singularity, the so-called "Big-Bang" singularity. Any hope of 
dealing with such singularities would be in the development of a concomitant and conducive quantum 
theory of gravity. In the absence of a full theory of quantum gravity, it would be useful to describe the 
quantum state of the universe within the context of quantum cosmology, introduced in the works of 
DeWitt [1]. In this formalism which is based on the canonical quantization procedure, the evolution 
of universe is described by a wave function in the mini-superspace. Two major approaches in this 
scenario are the tunneling proposal first developed by Vilenkin [2]- [6], and the no-boundary proposal 
of Hartle and Hawking [7]- [9]. In the tunneling proposal the wave function is so constructed as to 
create a universe emerging from nothing by a tunneling procedure through a potential barrier in the 
sense of usual quantum mechanics. In the no-boundary proposal of Hartle and hawking on the other 
hand, the wave function is constructed by a path integral over all compact Euclidean 4-manifolds. 
A problem related to this approach is that of signature transition from a Euclidean to a Lorentzian 
manifold. The notion of signature transition was first addressed in [7], where the authors argued 
that in quantum cosmology amplitudes for gravity should be expressed as the sum of all compact 
Riemannian manifolds whose boundaries are located at the signature changing hypersurface. In the 
traditional point of view, a feature in general relativity is that one usually fixes the signature of the 
space-time metric before trying to solve Einstein's field equations. However there is no a priori reason 
for doing so and it is well known that the field equations do not demand this property, that is, if 
one relaxes this condition one may find solutions to the field equations which when parameterized 
suitably, can either have Euclidean or Lorentzian signature [10]-[16]. 

An important ingredient in any model theory related to cosmology is the choice of the matter 
field used to couple with gravity. The most widely used matter source has traditionally been the 
perfect fluid. However, the scalar field has also been playing an increasingly important role in the 
more resent cosmological models. One main reason is that being a scalar field makes it somewhat 
easy to work with. Another matter field which has occasionally been studied in the literature is the 
massless or massive spinor field as the source of gravity. In general, theories studying spinor fields 
coupled to gravity result in Einstcin-Dirac systems which are not easy to solve. The cosmological 
solutions of such systems have been studied in a few cases by a number of authors [16]- [20]. The 
references [16] and [17] are notable in that the quantization of a spinor field coupled to gravity is 
studied in a Robertson- Walker background. 

In this paper we deal with the classical and quantum cosmology of a model in which a massive 
self interacting spinor field is coupled to gravity in a Bianchi type-I space-time. Bianchi models are 
the most popular anisotropic and homogeneous cosmologies and are studied in different levels by 
many authors, [19]- [26]. Prom the classical solutions of the resulting Einstein-Dirac system we have 
chosen those that exhibit a smooth transition from a Euclidean to a Lorentzian region, i.e., signature 
changing solutions. Although in classical gravity such solutions may not seem to be too interesting, 
we show that there is a close relationship between these solutions and the those resulting from the 
Wheeler-DeWitt (WD) equation in the corresponding quantum cosmology. It turns out that the WD 
equation possesses exact solutions in terms of confluent hypergeometric functions. We show that 
these solutions predict creation of a universe by a continuous transition from a classically forbidden 
(Euclidean) to a classically allowed (Lorentzian) domain and that they axe in agreement with the 
classical signature changing solutions. We also show that imposing boundary conditions on the wave 
function of the universe suggests a quantization condition for the mass of the spinor field. 



2 



2 Preliminary setup 



The Bianchi type I space-time is characterized in a comoving coordinate system by the metric 

ds'^ = -dt^ + a^{t)dx'^ + V + c^{t)dz'^, (1) 

where a(t), h{t) and c{t) are scale factors in x, y and z directions respectively. This metric is the 
simplest anisotropic and homogeneous cosmological model which, upon having equal scale factors 
becomes the flat Robertson- Walker metric. Such space-times have an Abelian symmetry group of 
translations with Killing vector fields ^ = {dx,dy,dz)- Of course, all the structure constants of such 
a symmetry group are zero. The scalar curvature corresponding to metric (1) is 

(a b c ah he ca\ 
a c ao be ca J 

where a dot represents differentiation with respect to t. We may parameterize the metric in such a 
way as to allow the Euclidean signature (+,+,+,+) becoming Lorcntzian (— , -|-, -|-, +) [11]-[15]. To 
this and other ends, we parameterize the metric as in [14] and [15] by adapting the chart {/?, x, y, z} 
where the hypersurface of signature change would be characterized by /3 = 0. The metric can then 
be parameterized in terms of the scale factors a(/3), b{(3), c(/3) and laps function (3 to take the form 

ds'^ = -I3d0^ + a^{j3)dx'^ + b^{l3)dy'^ + c^{f3)dz'^. (3) 

It is now clear that the sign of the laps function (3 determines the signature of the metric, being 
Lorentzian if /? > and Euclidean if < 0. For the Lorentzian region the traditional cosmic time 
can be recovered by the substitution t = Also, adapting the chart in this region 

we shall write any dynamical field such as the scale factors or matter fields as <I>(t) = <I>(/3(t)). As a 
result of the above discussion, we see that the signature changing hypersurface divides the manifold 
into two domains, Euclidean Ms and Lorentzian Mc- As has been discussed in [13], one may write 

Mer\Mc = ^, Mer\Mc = ^, Me^Mc = M, (4) 

where S represents the signature changing hypersurface, that is the hypersurface with /? = and M 
is the total manifold. From the point of view of the Einstein field equations, Euclidean and Lorentzian 
regions are classically forbidden and allowed solutions of the gravitational field equations respectively. 
We formulate our differential equations in a region which does not include /3 = and seek solutions 
for any dynamical field that smoothly passes through the /3 = hypersurface. These solutions are 
called signature changing solutions and we shall see that they are classical description of the quantum 
cosmological states of the model. Indeed, we will encounter the Euclidean and Lorentzian regions 
again when dealing with the solutions of the WD equation later on. 
To construct the field equations, let us start with the action 

{Lgrav + Lmatt)y/^d'^X, (5) 

where 

Lgrav =TZ- A, (6) 

is the Einstein-Hilbert Lagrangian for the gravitational field with cosmological constant A, and Lmatt 
represents the Lagrangian of the matter source, which we assume to be a massive spinor field. As 
is well known, the Dirac equation describing dynamics of a spinor field ^ can be obtained from the 
Lagrangian 



_ 1 

-'matt 2 



^j^^id^ + T^)^ - ^p0; - r^)7'^v] - v{^, V'), (7) 



where 7^ are the Dirac matrices associated with the space-time metric satisfying the Clifford algebra 
{t'^jT'^} = 2^'*^, are spin connections and V{^,il)) is a potential describing the interaction of the 
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spinor field with itself. The 7^* matrices are related to the flat Dirac matrices, 7", through the tetrads 



as follows 



(8) 



For metric (1) the tetrads can be easily obtained from their definition, that is, g^i, = e1^e^r]ab, leading 
to 

e« = diag(l, a, b, c), = diag(l, 1/a, 1/6, 1/c). (9) 
Also, the spin connections satisfy the relation 



r^ = ^ff.A(Va + raM<)7V- 



(10) 



Thus, for the line element (1), use of (8) and (9) yields 



To = 0, = - J7°V, = -^7V, Fa = —jV- 



(11) 



Here, 7° and 7* are the Dirac matrices in Minkowski space and we have adopted the following 
representation [27] 

f -i 0\ f a' \ 
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i 



7 



a' 



(12) 



The final remark about the Lagrangian (7) is that consistency of Einstein field equations with a spinor 
field as the matter source in the background metric (1) requires the spinor field to be dependent 
on t only, that is i/j = tp{t) [16]. 

The preliminary set-up for writing the action is now complete. By substituting (2), (6) and 
(7) into (5) and integrating over spatial dimensions, we are led to an efi^ective Lagrangian in the 
mini-superspace {a, 6, c, V', V'} 



£. = dbc + abc + dbc + Aabc + -a6c[^7°V' - "07° - 2F(V', V")]- 



(13) 



3 Field equations 

Variation of Lagrangian (13) with respect to ■0, ijj, a, b and c yields the equations of motion for the 
spinor and gravitational fields respectively, that is 



• 1 (a b c\ , ^ 



J fa b c\ - dV Q 



b c be ^ 1 
6 + c + 6^-^=2 



c d ac , 1 

-H \ A = - 

c a ac 2 



d b db , 1 
a b ab 2 



-dV ^ dV^ 

dtjj dip 

'-dV ^ dV ^ 
dip dtp 

'-dV ^ dV ^ 
dip dip 



0, 

-ViiP,iP), 

-v{iPA), 

-V{iP,iP). 



Also, we have the "zero- energy" condition given by 



dC. dC- dC. dC ■ jdC , ^ 

W = + + + -jiP + V— - £ = 0, 
da db dc dip dip 



(14) 

(15) 

(16) 
(17) 

(18) 
(19) 
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which yields the constraint equation 

ab be dc . ^^,t ,^ 

- + - + --A = -y • 20 

ab be ae 

It is clear that the right hand side of equations (16)-(18) and (20) represents the components of the 
energy-momentum tensor where the energy density of the spinor field is given by 

/9 = roo = -y(^,V). (21) 

Finally we write the Hamiltonian constraint (19) in terms of the momenta conjugate to our dynamical 
variables, giving 

5£ ; , . dC . . dC 1 . , 

Pa = -^ = oc + bc, pb = —f = ac + ac, p^ = — =ah + ab, (22) 
oa db oc 

and also ^ ^ 

p^ = —j = -abei!-f^, pT = — = --abej^^. (23) 
d^p ^ (9V 2 

In terms of conjugate momenta, the Hamiltonian is given by 

l/^a2,^2, C2\l fPaPb , PbPc . PaPc\ k , . , MfT l\ n fOA\ 

In the next section we shall present the solutions of the classical field equations (14)-(18). Although, 
these equations can be solved after a suitable form for the potential V{'4),ij)) has been chosen, see 
[19] and [20], Hamiltonian (24) has not the desired form for the construction of the WD equation 
describing the relevant quantum cosmology. Thus, to simplify Lagrangian (13), consider the following 
change of variables [28] 

= 6"+"+^"', 6 = 6"+''-^'", c = e"-2^ (25) 
In terms of these variables, Lagrangian (13) takes the form 

£ = 3 (li^ - - v?^ e^" + Ae^" + ie^" [^^7°^' - V'tV - 2t^(Vi, V-)] , (26) 
with the corresponding "zero-energy" condition 

e-3«7^ = 3 (^^(2 _ ^2 _ ^2^ _ ^ ^ yl^^^ ^) ^ Q_ (27) 

Also, the momenta conjugate to u, v and w are 

BC BC BC 

Pu = ^ = Que''', Pv = ^ = -Gve'\ P-^ = ^ = "^^e^", (28) 
ou ov ow 

which give rise to the following Hamiltonian for our dynamical system 

n = ^e-3« [pI -pI- pI) + [V{^, i;) - A] e^" = 0. (29) 

Now, variation of Lagrangian (26) with respect to its dynamical variables yields the following field 
equations 

^ + ^^^ + /|K=o, (30) 
^ + 1^7° = 0, (31) 

2u + 3u^ + 3{v^ + w'^)-A-^ [v^tV - V'tV - 2F('iA, tp)] = 0, (32) 

(i,e3"y = 0, (33) 

(u;e^")=0. (34) 

Our goal would now be to concentrate on the solutions of these equations for a suitable form of the 
potential V. 
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4 Classical solutions 



The system of differential equations (30)- (34) are the Einstein-Dirac system for the gravitational field 
coupled to a massive self interacting spinor field. Since integrability of this system directly depends 
on the choice of a form for y(^,'i/'), it is appropriate to concentrate on this point first. Interesting 
forms for such potentials should involve terms which would describe self-interacting spinor fields. 
The potential is usually an invariant function constructed from the spinor ip and its adjoint ip and 
we require it to satisfy the following property 

This property is a mild restriction on the potential as may easily be verified for the most commonly 
used ones. One such form for the potential that has property (35) is 

F(^, = m^pi; + A(V^^/')^ (36) 

where m is the mass of the spinor field and A is a coupling constant. With this potential we can 
immediately integrate equations (30) and (31) to obtain -0^ = Ce~^", where C is an integrating 
constant. Bearing in the mind from equation (21) that the energy density of the spinor field should 
have a positive value and with the use of (36) we are led to a negative value for C which we take 
choose as —1, hence 

= (37) 
The next step is to solve equations (27) and (32)- (34). With potential (36) these equations read 

V = cie"^", (38) 

w = C2e-^", (39) 

2u + Sii^ + — e-'5" - A = 0, (40) 

O 

3u^ - - A - = 0, (41) 

where in deriving them we have integrated equations (33) and (34) and substituted the results in (27) 
and (32). ci and C2 are two integrating constants and = 9(cf + C2 — A/3). Of course, equations (40) 
and (41) are not independent and a solution to one should satisfy the other automatically. Integration 
of equation (41) would be made easier after a change of variable r = e^", resulting in 

= SAt^ + 3mr + C^. (42) 

This would now be a simple equation to integrate. Depending on the sign of A, there are three classes 
of solutions for A = 0, A < and A > 0. We are interested on those values which would allow for 
signature changing solutions to be constructed. For A = 0, we have 

3m / ,s9 / 

where ^ is a constant. It is clear that this solution, after substitution t = does not exhibit 

smooth transitions from /3 < to /? > regions. In the case of A < 0, integrating equation (42) 
results in 



m / m 



2 



C2 
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where the integrating constant is chosen so that f(0) = 0. The functions u, v and w can then be 
easily obtained in terms of the evolution parameter /? as 



u{(3) = In 



m 



+ 



C2 



4A2 -3A 



cos 



1/3 



w{P) = I In 



{B - A) tan (^^^A/SP^/^^ + C/V^3A 



{B - A) tan (^^^A/SP^/^) - C/V^SA 
(B - A) tan (y^A/S/J^/^) + C/V^ 



{B - A) tan (y^A/S/J^/s) _ C/V^ 



(45) 

(46) 
(47) 



where A 



and B 



4A- 



+ 



C2 



with the scale factors being obtained from (25). A glance at 
equations (45)- (47) shows their regular behavior both in Euclidean and Lorentzian regions passing 
smoothly through signature changing hypersurface /3 = 0. We should note that to avoid imaginary 
values for n(/3) in (45), we must restrict the allowed values of (3 to the interval < /3 < /3o = (rfx)^^"^- 
However, this is not a sever restriction since /? = is the beginning of the Lorentzian region and with 
the present bound on the cosmological constant |A| ~ 10~^^ cm~^, one has Pq >present age of the 
universe. In the limits of /? — ^ and /3 — Pq, equations (46) and (47) show that v{f3),w{f3) — 0, 
i.e. the scale factors a, b and c become equal and the universe approaches that of the flat Robertson- 
Walker in these limits. In the case of A > the solutions can easily be obtained by the replacement of 
the "cos" function in (44) with its hyperbolic counterpart and a quick look at the resulting solutions 
shows that they do not exhibit the signature changing behavior, see [16] for details. 

Thus, in summary, the above discussion shows that within the context of our model, a universe 
with negative cosmological constant can undergo signature transition from a Euclidean to a Lorentzian 
domain through the /3 = hypersurface, where its matter source is a spinor field. Indeed solutions 
(45)- (47) describe the classical evolution of a universe from a classically forbidden region where there 
is no time, to a classically allowed Lorentzian region. This issue will become more clear when wc 
study the quantum cosmology of the model in the next section. Figure 1 shows the behavior of the 
scale factors a(/3), 6(/3) and c(/3) as (3 changes from negative to positive values, showing a change of 
signature. 



5 Quantum cosmology 

The study of quantum cosmology of the model presented above is the goal we shall pursue in this 
section. For this purpose we quantize the dynamical variables of the model with the use of the WD 
equation, that is, Ji^ = 0, where Ti is the operator form of the Hamiltonian given by equation (29) 
and * is the wave function of the universe, a function of the scale factors and matter fields. With the 
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replacement Pu —>■ and similarly for py and py,, and choosing potential (36), the WD equation 

reads 



^ ^ ^] 12 r 

du^ dv"^ dw'^ I ^ 



-me 



3n 



Ae^" + A 



'^{u,v,w) = 0, 



(48) 



where we have used equation (37) to reduce the degrees of freedom in the wave function. The solutions 
of the above differential equation are separable in the form ^i{u,v,w) = U {u)V {v)W {w) , leading to 



1 (PV 

V dv^ ' 

1 (fW 
W~duP 



+ 12 (-me^" - Ae^" + A ± .r^) C/ = 0, 



(49) 

(50) 
(51) 



where rf and Q"^ are separation constants and 



12 



Before presenting the solutions of the 



above equations, note that equation (51) is a Schrodinger-like equation for a particle with zero energy 
moving in the field of the potential 



U{u) 



-Ae 



Gu 



me 



3m 



+ A±?^. 



(52) 



For a negative cosmological constant this potential has a minimum 

at = ln(^)V3. In the presence 
of this potential the mini-superspace can be divided into two regions, U > and U < which could 
be termed as the classically forbidden or Euclidean and classically allowed or Lorentzian regions 
respectively. The boundary between the two regions is given by W = 0, that is at 



u± 



In 



m 



4A2 



-A 



1/3 



(53) 



Comparison of the above values of u and the classical solutions (45) suggests that they correspond 
to /3 = and (3 = [3q. Thus, the same boundary separates the Euclidean and Lorentzian regions 
in both classical and quantum solutions. In the Euclidean domain we have the wave-functions with 
exponential behavior and in the Lorentzean region we have the wave-functions of oscillatory nature 
[3, 4]. Prom equations (45)-(47) it is clear that we have no classical solutions that extend to infinite 
values of u, v and w. Also the beginning of evolution in the Lorentzian region is represented by 
V = w = Thus we are led to the following boundary conditions 



^' 0, as u,v,w ^ +00, 
^' — > 0, as v, 1/; — — GO. 



(54) 



The above discussion suggests that we take the right hand side of equations (49) and (50) with 
the upper sign, leading to solutions 



in the Euclidean region, and 



V{v) 



V{v) 



-\r)\v 



W{w) 



W{w] 



(55) 



(56) 



in the Lorentzian region. Now let us deal with the solution of equation (51). This equation, after a 
change of variable r = cje^" with u) = 4i/— A/3 and transformation U = {t/u)'^^'^^, becomes 



+ 



-1 K 1/4 

-r + - + 



0, 



(57) 
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where k = m/\/— 3A and fi^ = |(A + i;^). The above equation is the Whittaker differential equa- 
tion and its solutions can be written in terms of confluent hypergeometric functions M{a, b; x) and 
U{a,b;x) as 

0(r) = [cUip - K + 1/2, 2^ + 1; r) + c'M(^ - k + 1/2, 2^ + 1; r)] . (58) 

Since the asymptotic behavior of M{a,b;x) ^ e^/x*~" [29], we take d = 0. Therefore 



A feature of the quantum cosmology of our model is the quantization condition for the mass of the 
spinor field. A glance at equations (57) shows that its solutions have a behavior of the form (f> ~ e""^/^ 
when T +00 (or u +oo) and (j) ^ m the limit r ^ (or u — oo). Demanding the same 
limiting behavior for solution (59), one can easily see that the function U{a,b;x) should reduce to a 
polynomial and this happens when a = —n [29]. This yields a quantization condition for the mass of 
the spinor field as follows 

m = ^ ^2n + 1 + 4y^ ^"^^^ j w. (60) 
Thus, the above discussions lead us to the following eigenfunctions 



(61) 



with the upper and lower signs signifying the Euclidean and Lorentzian regions respectively. The 
general solution of the WD equation can then written as 

■^{U,V,W) =^Cn^n{u,V,w). (62) 
n 

To have an exponential wave function in the Euclidean domain wc must take n = in (62), i.e. the 
matter must be in its ground state. Taking more terms in (62) leads us to the wave function in the 
Lorentzian region. Summarizing, we have the following wave functions 



(63) 



in the Euclidean region and 



(64) 



in the Lorentzian region. The wave functions (63) and (64) describe a universe emerging out of 
the Euclidean region with a smoothly changing signature and correspond to the signature changing 
classical solutions. The creation of the Lorentzian universe in this scenario which is characterized 
by the smooth passage through the Euclidean region is comparable to the quantum tunneling from 
nothing in the Vilenkin's proposal [2]-[6] where nothing is a 3-geometry of vanishing size or a point. 
However, one should note that potential (52) has no maximum and therefore should not be considered 
as a potential barrier like those described in [2]- [6]. Since the potential (52) has a dip the wave- function 
(64) in the Lorentzian region resembles bound states of the Schrodinger equation in a potential well. 

In the case of a positive cosmological constant, potential (52) is a monotonically decreasing func- 
tion of u. Although in this case the mini-superspace can be divided into two regions characterized 
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hy U > and U < too, in the Lorentzian region, U < 0, we have solutions which extend to infinite 
values of u, instead of being oscillatory as in the case of A < 0. Thus, there is no mechanism for a 
smooth transition from the classically forbidden to the classically allowed regions and the creation 
of a Lorentzian universe would not be possible. We have seen before in section 4 that the classical 
solutions with A > do not exhibit any signature changing properties either, that is, the classical 
and quantum solutions are in agreement in this case as well. 



6 Conclusions 

In this paper we have studied the classical and quantum evolution of the cosmological solutions of the 
Einstein-Dirac system in a Bianchi type I background. Our approach is a dynamical system approach 
in which the mini-superspace is constructed from three scale factors and components of the spinor 
field as the matter source. From the classical solutions of this system we have chosen those that admit 
a degenerate metric for which the scale factors of the universe have a continuous behavior in passing 
from a classically forbidden (Euclidean) to a classically allowed (Lorentzian) region. We have shown 
that this happens when the cosmological constant is negative. The corresponding cosmology begins in 
the Lorentzian domain by admitting a flat Robertson- Walker metric, evolves according to equations 
(45)- (47) and finally approaches a flat Robertson- Walker universe again. The quantum cosmology 
of the model presented above and the ensuing WD equation is amenable to exact solutions in terms 
of confluent hypergeometric functions. We have found that the division of the mini-superspace into 
the classically forbidden and allowed regions also happens at the quantum level. In the case of a 
negative cosmological constant the behavior of the wave-functions are exponential in the Euclidean 
domain and oscillatory in the Lorentzian region. These wave-functions correspond to the classical 
signature changing solutions and as such, could be useful in understanding the initial conditions of 
the universe. Within the context of this model the creation of the universe is described by a smooth 
transition from a Euclidean to a Lorentzian region. Another feature of the quantum cosmology of 
our model is a quantization condition which leads to a spectrum for the allowed values of the mass 
of the spinor field. 
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